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O We adopt a t\ — 12 — t^ — J — G model for explanation of x = 1/8 anomaly in 
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La2_ I Sr :r Cu04 family compound. The calculated charge susceptibility shows 
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a maximum near (n,Tr) at intermediate temperatures and near (tt,tt/2) as 
temperature approaches zero, in agreement with neutron scattering experi- 



o 
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ments. Coulomb repulsion G between the first neighbors turns out to be the 
source of Charge Density Waves (CDW) in narrow band ,t%~ < G. 
. For physically realistic hopping values we obtain the CDW amplitude eQ = x. 
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The in-phase domain structure as a candidate for "stripe" picture is proposed. 
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It is widely accepted, that high-temperature superconductivity (HTSC) appears in the 
vicinity of metal- insulator transition (MIT) in layered copper oxides |I| . In addition, these 
compounds possess a pseudogap phenomenon in the normal state of underdoped regime, 
namely the density of states suppression at Fermi level. The superconducting CuC>2 planes 
in La 2 _ a; Sr a .Cu04, YBa2Cu 3 07_ y and other copper oxide layered HTSC compounds present 
nearly an ideal realization of two-dimensional (2D) systems. Low dimensional character of 
these compounds suggests, that various phases are possible due to instabilities. It is known, 
that s-wave BCS Superconductivity compete with Charge Density Waves (CDW), where 
the anomalous averages were due to Frohlich interaction and did not depend on the wave 
vector, ci-wave superconductivity was also found to compete with id-wave Peierls instability 
in bilayered underdoped HTSC, where both gaps were due to Superexhange and Coulomb 
interactions between the nearest neighbors ||. Very recently this sliding vortex-antivortex 
current structures in Cu0 2 planes attracted a renewed interest in connection with unusual 
nonsuperconducting state properties of layered copper oxide HTSC compounds @-§|]. 

Strong depression of T c in La2_ x Ba x Cu04 [J5], 

La2_ x Sr x Cu0 4 |J and Lai. 48 Ndo.4Sr .i2Cu04 ]7] within narrow region of x at x ~ 1/8 is 
the well known important problem of the doped Cu0 2 plane in the HTSC. The discovery of 
stripe like structure via neutron diffraction in Lai.4sNdo.4Sro.i2Cu04 |§ gave a new acceler- 
ating impulse for study this intriguing problem (see |§ for review). New ideas regarding the 
relation between the fluctuating stripe motion and HTSC phenomenon were speculated |10 



The Q dependence of peaks as seen in neutron diffraction experiments seems to provide an 
evidence that stripes are a kind of CDW instability; this scenario was considered in a number 



of papers (see e.g. JTl|). However, the crucial question, namely why at so narrow interval of 
concentration around x ~ 1/8 the mystery of Tc strong depression occurs and accompanies 
by stripes pattern phenomenon is not been completely understood. Recently this problem 



was investigated by magnetic resonance method |12| and new evidence of the quasi-static 



charge-spins ordering at x ~ 1/8 has been found at low temperatures in La2_ x Ba x Cu04 [T3 



and in Lai.4sEuo.4Sro.i2Cu04 |141| . Note that x-ray diffraction measurements have confirmed 



the existence of charge order in Lai.4sNdo.4Sro.i2Cu04 |L5] . 



In this Letter we will concentrate ourselves to the description of properties for particular 
underdoped case x — 1/8 in La2- x Sr~Cu04 on the microscopical basis. We start from the 
t\ — ti — ts — J — G model Hamiltonian 
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where 5i = 1 — rii is the number of extra holes per unit cell, iJrf'W'*) are creation 
(annihilation) operators constructed on the singlet combination of copper (d) and oxygen 
(p) states basis (see Ref. |16| for details). The last term takes into account the polarization 
of nearest neighbor (n.n.) copper spins (p runs over n.n.) around copper-oxygen singlet at 
site i as described in |I7| . 

Following the idea of Zhang and Rice |18|] about copper-oxygen singlets formation it was 



shown ||T6||, that it is possible to describe correctly the elementary excitations spectrum in 
cuprates. This singlet correlated band is analogous to upper Hubbard band with essential 
distinction - the subband splitting is much smaller, compared to Hubbard model. Therefore 
it is possible to apply Hubbard formalism without strict restriction on t and J values in 
t — J model. According to Angle Resolved Photoemission Electron Spectroscopy (ARPES) 



the bandwidth W ~ 3 J for optimally doped HTSC [19]. In undoped insulating Sr 2 Cu0 2 Cl2, 



W ~ J, in agreement with Monte-Carlo and Exact Diagonalization simulations |19| , and 
the valence-band dispersion has a maximum located near (7r/2,7r/2) f20fl . The value of bare 
hopping between the first neighbors ti = 78 meV is known from comparison with ARPES 
P~6|1 . We stress the importance of t 2 and £3 hopping integrals, that expand our model and 
strongly affect on the Fermi surface form. The value of superexchange interaction between 
the nearest Cu neighbors for La2_x.Sr2.CuO4 is known from neutron scattering data J=135 



meV P_|| . These arguments enforces us to study the opposite region of values (ti, t 2 ,t 3 < J), 
compared to usually used in celebrated t — J model. 

The hoppings Uj are strongly affected by electron and especially by spin-spin correlations, 



resulting to effective values of hoppings 

feff -l p 
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for n=\j — z|-th neighbor. P p d=(l + x) /2 is Hubbard electronic reduction factor. The 



spin-spin correlation functions (SjSj) have been calculated in a variety of approaches [22 



and numerical methods [T]J within the 2D t — J model. The most notable result is that 
(SjSj)i ~ —0.333 in undoped antiferromagnet at zero temperature. These calculations 
suggest also that with doping the spin-spin correlations goes down first slowly in lightly 
doped and underdoped HTSC oxides and then rapidly decreases for optimally doped and 
overdoped regime. Therefore, for i = 1/8 we set (SjSj)i = —Pp d ~ —0.3164 so that t± = 0. 
Strong antiferromagnetic correlations block the hopping between the first neighbors so that 
the order is undamped H. Our initial Fermi surface with dispersion 



4t%P cos k x cos k v + 2t^f {cos 2k x + cos 2k 



(3) 



has pockets centered around (7r/2,7r/2) and is similar to that observed by LaRosa et al. [20 
in S^CuC^Ck). 

To select the instability wavevector for charge subsystem we have calculated the charge 
susceptibility, using the expressions given in P3[ . The results are shown on Fig. 1. It is seen, 
that at low temperatures charge susceptibility has a maximum near (vr/2, tt), in remarkable 
agreement with neutron scattering experiments || . With increasing temperature the charge 
susceptibility shows a peak around Q = (7r,7r). These results are stable with respect to 
moderate variation of G value. 

To formulate the self-consistent system of equations we use Green's function method 



and apply Roth's decoupling scheme ||24|| . The elementary excitations spectrum of the 
Hamiltonian ([]]) is given by the equation 
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The quasiparticle dispersion relation has the form 
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where £ k = 2ii (cos k x + cos fc^) + 4i 2 cos /c x cos hy ~t~ (cos 2/u^ ~\~ cos 2/cy ) . Tlie order parameter 
q is responsible for abnormal properties of HTSC layered copper oxides due to instabilities. 
It is caused by short ranged interactions and given by 
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where J k = 2 J (cos /c x + cos k y ), G k — 2G (cos fc x + cos k y ) and a = ± is spin index. The 
latter contribution in Eq. (|j), GqCq, will be important and we note, that Gq/G = —4 is 
negative due to instability vector Q = (tt,7t). The eigenvalues of the matrix (|J) are given 
by 



E 



J k,2k _ 2 V^k+Q + £ kj A 1 =t 2 -^12 



(7) 



where 



#f 2 = V £ k+Q - ^k + M+Q.Q^Q- 



(8) 



It is seen that in the = case r/ kj Q = ?7 k+ Q_Q an d the quasiparticle energy is real, i.e. 
the regime realized has no damping. 
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are Fourier amplitudes of CDW and Spin Density Wave (SDW) respectively and can be 
expressed via thermodynamic averages as 

1 



(eq) 
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We note, that the id-wave order does not contribute to expectation values of CDW (ecj) and 

sdw 

The relevant thermal average have the following form: 



x-2cr vcr2\ _ "pd^k.Q 



L k+Q vi k / rpa 
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In the ti = 0, t 2 ,£3 -C G limit suggestions on solution of system of equations (|q-|12|) 
may be obtained even analytically. Suppose, that we have non-negligible CDW amplitude cq. 
Then, due to large G, compared to hoppings, value, the contribution GqCq will dominate 
in (§). Next, we note, that the rj order parameter will dominate in Eq. @ and hence 
E° 2 « |2^q|. Therefore the absolute value of thermal average ( |T2"D is equal to |-P p <f- 



Substituting this result to Eqs. ([10], |TT]) we obtain, that the SDW amplitude is sq = and 
CDW amplitude cq = \P p d- 

Despite many attempts to calculate the numerical values of hoppings t^f and t^f , their 
values are precisely still unknown together with the spin-spin correlation functions between 
the next-next nearest neighbors. The only known result from ARPES experiments and exact 
numerical simulations within 2D t — J model is that the bandwidth W ~ J for undoped 
antiferromagnet, that gives a restriction to the t% and t$ values. We set the values of 
effective hoppings \t% \ = 10 meV and £3 = 20 meV, that are, in our opinion, the most rele- 
vant values for La2_ 2; Sr a ;Cu04 with x — 1/8. For the spin-spin correlation function between 
the second neighbors we set (SjS 3 -)2 = 0.15, due to reasons described above. The value 
and even the sign of next-next nearest neighbors spin-spin correlation function is completely 



unknown, therefore we set (SjSj)3 = 0. F{ ~ 4 meV was estimated in 17]. The value of 



screened Coulomb repulsion at x ~ 1/8 is G ~ 75 meV p3 |. We have solved the system of 
equations (|6|-|P2|) with this set of parameters. We have obtained, that in the low temperature 
limit the CDW amplitude eg = x, where x is the number of holes in CUO2 plane due to Sr 
doping. In this case the undamped CDW completely destroy superconductivity at x — 1/8 
in La2-a;Sr a ;Cu04 and are responsible for drastic change in temperature dependencies of 
various physical quantities. The calculated wavevector dependence of 7/£ order parameter at 



T= 50 K is presented in Fig. |2|. It is seen, that its value is mainly determined by wavevector 
independent CDW contribution GqCq. The temperature dependence of CDW amplitude eg 
determined by temperature dependence of thermal average fll2"|) , which are, in turn, deter- 
mined by mean field approximation via Fermi distribution functions f (E^) , f (E^) . The 
calculated CDW closure temperature is Tcdw — 350 K. 

The described above scenario explains the reason of anomaly at x = 1/8. At this doping 
level conventional CDW are allowed (without damping) and they suppress superconductiv- 
ity. According to our calculation CDW regime appeared at x = 1/8 far above To The 
experimental observation of critical temperature Td2 — 138 K in La2_ x _ 3/ Eu y Sr :r Cu04 with 



x = 1/8 at which the copper relaxation rate 1/Ti suddenly goes down [|13j can be naturally 
interpreted as a transition related to CDW state. With lowering temperature the pinning 
process due to freezing of the spin fluctuation becomes important and charge-spin ordering 
tends to be quasi-static. In Fig. |3] we show our variant of copper - oxygen singles ordering 
in plane which we expect to be appeared in La2_ a; _j / Rej / Sr a ,Cu04 (Re=Nd and Eu) below 
T ~ 9 K. The displayed stripe pattern has many common features with that suggested by 
Tranquada et al. ||. The distinction is that our picture corresponds to the so-called in- 
phase domain walls in contrast to antiphase domain proposed in ||. As one can see our 
model explains the observed magnetic superlattice peaks of the type (1/2, ±e), (1/2,0) (in 2ti 
units) and charge-ordered peaks at (2e,0,0) |§, too. In addition it is naturally follows from 



widely accepted copper-oxygen singlet formation picture |L8|], which provides the minimum 
in the energy of the exchange interactions between copper - oxygen and copper-copper spins 
||17|| . The last term in ([I]) plays crucial role in this aspect. It was pointed in Ref. that 



in-phase domain even better reproduce the observed high frequency tail in NQR spectra in 
La2_ x Ba a: Cu04 with x = 1/8 at low temperatures compared to anti-phase domain model ||. 
This conclusion is consistent with the ferromagnetic stripe ordering suggestion |25] based 
on strong sensitivity to magnetic fields. The Coulomb term (G) is important, since without 
it there is no CDW-like instability. This statement is in agreement with conclusion of Eft 
that pure t — J model does not give any stripes. 



Finally we comment the results of huge isotope shift changes with Sr doping in 



La 2 _ x Sr :r Cu04 According to our findings the competition between ci-wave superconduc- 
tivity and CDW is stronger, compared to competition between (i-wave SC and order, caused 
by short range interactions. Moreover, the latter type of competition becomes stronger with 
decreasing of doping. Thus we are able to explain the gradual increase of isotope effect with 
decreasing of Sr doping in underdoped regime together with nearly divergent behavior of 
isotope exponent near x = 1/8 anomaly in L^-xSr^CuC^. 

In summary, the calculated charge susceptibility shows a maximum near (tt, it) at in- 
termediate temperatures and near (71,11/2) as temperature approaches zero, in agreement 
with neutron scattering experiments. We have shown, that in two-dimensional systems with 
narrow bands, so that the hoppings t\ , < G, t± = due to strong antiferromagnetic 
correlations, and with instability antiferromagnetic wave vector Q = (ir,ir), due to large 
Coulomb repulsion between the first neighbors G the CDW are formed. For physically real- 
istic values tif ,£3 at x = 1/8 in La2-xSr x Cu04 the CDW amplitude cq = x, where x is Sr 
content. Based on the calculated charge susceptibility peaks position the in-phase domain 
structure as a candidate for "stripe" picture is proposed. In our opinion, this scenario sheds 
light on drastic change of various physical quantities at "magic" value of Sr content x = 1/8 
in La 2 _ a; Sr a .Cu04. 
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FIGURES 




FIG. 1. The calculated real part of charge susceptibility at T = 10 K with the parameters as 
described in the text. 




FIG. 2. The calculated wavevector dependence of rfa order parameter. The main value of t/£ 
originates by wavevector independent Coulomb x CDW amplitude contribution Gqcq. 
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FIG. 3. Sketch of in-phase "stripe" picture. The circles correspond to oxygen holes. 
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